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LONG-TERM PREDICTIVE CAPABILITY OF EROSION MODELS 
P. Veerabhadra Rao and Donald H. Buckley 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

ABSTRACT 

This paper reports a brief overview of long-term cavitation and liquid 
impingement erosion and modeling methods proposed by different investigators, 
incliding the curve-fit approach recently suggested from this laboratory. A 
table is prepared to highlight the number of variables necessary for each model 
in order to compute the erosion-versus-time curves. A power law relation based 
on the average erosion rate is suggested which may solve several modeling 
problems. 


INTRODUCTION 


Long-term prediction of erosion due to cavitation and liquid impingement 
has become very important in view of severe erosion problems associated with 
hydraulic turbines and pumps due to cavitation, with aircraft surfaces due to 
rain drops, and with steam turbine blades due to impingement of condensed drop- 
lets. Extended periods of reliable operation in all situations reemphasized 
the necessity for erosion-free performance or the alleviation of erosion. In 
most cases complete elimination of erosion is not possible. Hence, it is 
necessary to establish the total erosion of a material for prolonged operation 
so that the particular component may be changed, or so that highly resistant 
materials may be used to increase component life. 

Although Honegger (ref. 1) was the first investigator to notice the effect 
of exposure time on erosion rate, Fyall el al . (ref. 2) in 1957, Hobbs (ref. 3) 
in 1962, and Thiruvengadam (refs. 4 and 5"]" Tn the early I960' s clearly observed 
the influence of time on instantaneous erosion rate. Investigations by 
Thiruvengadam and Presier (ref. 6), Plesset and Devine (ref. 7), Heymann 
(ref. 8), and Tichler and de Gee (ref. 9) have become classic studies. There 
have been, however, several discrepancies in the agreement of the type of 
erosion-rate-versus-time curves. 2 Using shapes of curves obtained earlier, 
many models and formulations have been presented by different investigators 
for the long-term prediction of cavitation and liquid impingement erosion 
(refs. 8 to 14). The details of the models are reviewed and presented in 
reference 15. Table 1 presents models proposed and the number of variables 
necessary to predict erosion using these models (refs. 14 and 15). 

# 


^Instantaneous erosion rate equals slope of local tangent on erosion-versus- 
time curves 

2 Curves reported in reference 6 contain incubation, acceleration, decelera- 
tion, and steady-state zones; in reference 7, incubation, acceleration, steady- 
state, and deceleration zones; in reference 8, peak erosion and deceleration 
zones or acceleration zone and several cycles of peaks of erosion rate; and in 
reference 9, incubation, acceleration, first steady-state, deceleration and 
second steady-state periods. 
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A constant or coefficient 

ai parameter representing inner friction of material 

during plastic deformation 

a 2 a scale of cavitation strength of material 

B exponential constant, (fig. 1 ) or 

rcoeff icient (eq. (4)) 

C coefficient (eq. (7)) 

e erosion rate at any time x 

• 

ej„ ax maximum or peak erosion rate 


T 



max 


I e intensity of erosion at any time x 

Wx maximum intensity of erosion 


F (dn/dx) , 

x® i 

max maximum 

N e generalized nondimensional erosion resistance 

parameter defined as the ratio of the measured 
erosion rate of the material to the measured 
erosion rate of the material to be evaluated 
n attenuation exponent (fig. 2) or exponent 

(eqs. (1) and (2)) 

n w number of craters per unit area in final steady- 

state period 
P exponent 

R averaged erosion rate up to a cumulative material 

loss Y 

R c resistance against cavitation erosion under the 

hydrodynamic conditions as in magnetostriction 
oscillator 

Ry averaged erosion rate at tangent point corre- 

sponding to a cumulative material loss Yj 
R lu resistance against cavitation erosion during the 

second steady-state condition 
r mean depth of erosion at any time t 

r^ mean depth of erosion at which the effect of 

crater formation becomes manifest 
t exposure time to cavitation (figs. 2, 4, and 5) or 

incubation period (fig. 3) 

t, n time corresponding to the maximum intensity of 

erosion 

t 0 incubation period (fig. 4) 

tw time required to reach a mean erosion depth of 

Y (sec) 
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U a volume of liquid impinging per unit area per 

unit time (pm/sec) 

v cavitation damage rate - Am/ a t 

V cumulation volume loss, mm j 

V 0 normal component of impact velocity (m/sec) 

Vg steady-state cavitation damage rate 

Y cumulative mean depth of erosion or material loss 

at any instant 

Yy cumulative mean depth of erosion or material loss 

at tangent point 

o Wei bull shape parameter (fig. 2) or parameter 

representing inner friction of material during 
plastic deformation « a^/2 (fig. 5). 

B parameter representing square root of a scale of 

cavitation strength of material 

am cavitation damage of a material during time at 

at time increment 

6 parameter expressing cavitation properties of the 

material = o/e 

e a normalized average erosion rate with respect to 

peak erosion rate 


n efficiency of erosion at any time t * 1 - exp (-*<*) 

^ n ' n max 

n max efficiency of erosion corresponding to the peak 

or maximum intensity I max 

0 limiting value of 5 

k a proportionality constant symbolizing the 

increase in mean depth of erosion (r _> rb) 
which would be necessary to form n°° craters 
v relative cavitation damage =* Am/(v s At) = v/v s 

1 log t 

t t/t^ (in Thiruvengadam's theory of erosion, 

fig. 2) or exposure time (according to Hoff and 
Langbein equation, fig. 3) or relative 
cavitation time (fig. 5) = 0 t 
ratio between the rate of erosion in the final 
steady-state period and the rate of erosion in 

the first steady-state period 


LONG-TERM PREDICTION MODELS 


Equations or models for the prediction of erosion rate with respect 
to time have been proposed by Heymann (refs 8 and 12), Tichler and de Gee 
(ref. 9), Thiruvendagam (ref. 10), Hoff and Langbein (ref. 11), Noskievic 
(ref. 13), Engel (ref. 16), and Perelman and Denisov (ref. 17). McGuiness and 
Thiruvengadam (ref. 18) have studied the influence of corrosion on erosion in 
their modeling efforts. Others (refs. 19 and 20) have presented nomograms and 
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graphs. A brief description of the important contributions of these models 
is outlined below in order to explain the current status of erosion-rate- 
versus-time predictions for long-term exposures. 


Heymann's Models 


Elementary model . - Heymann (ref. 8) developed an elementary statistical 
erosion-rate-versus-time model for liquid impingement and cavitation erosion 
conditions of different materials wherein fatigue is the predominant failure 
mechanism. The model requires four parameters to obtain instantaneous erosion- 
rate-versus-time curves (table I). The fit of experimental data in certain 
real situations is very convincing with the use of normal distributions trun- 
cated and normalized over a finite time span. 


Elaborate model . - Heymann's elaborated model (ref. 8) permits the speci- 
f i cat ion of a different distribution function for each level below the original 
surface and of two different functions for the original surface. In this model 
the log-normal distribution is adopted. The inclusion of the median lifetime 
for the unaffected surface has significantly improved the predictions. 


Curve-fit approach . - Heymann (ref. 12) suggested a simple curve— fit approach 
"using tangent (cumulative average) rate of erosion 3 and tangent mean cumula- 
tive depth of erosion (volume Toss) to predict the erosion rate which follows 
the peak erosion rate (fig. 1). Equations suggested for the calculation of 
normalized average erosion rate (R/Rj) and for time (t y ) to reach a mean 
erosion rate require three and four parameters, respectively, to compute the 
erosion-versus-time history (table I). For particular liquid impingement and 
cavitation erosion data sets, this approach appeared promising (ref. 12). 


Thiruvengadam' s Nomogram and Theory of F.rosion 


Strain energy . - Using a strain energy theory, 4 a nomogram was developed 
based on cavitation erosion data (refs. 19 and 21). This nomogram has been 
used by design engineers to predict life of materials with a knowledge of 
erosion intensity and strain energy. Unfortunately, strain energy is a good 
predictor only for highly ductile materials. 


Concept of erosion strength and theory of erosion . - In view of the limita- 
tions of strain energy, Thiruvengadam later developed the concept of erosion 
strength 5 . Using this concept a theory of erosion (ref. 10) was developed 


^Cumulative average erosion rate equals slope of line joining origin and 
point of consideration on erosion-versus-time curve. 

4 The area of a stress-strain curve is a measure of this energy per unit 
volume. 

5 The energy-absorbing capacity of the material per unit volume under the 
action of erosive forces or the ratio of energy absorbed by the material 
eroded to the volume of material eroded. 
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to predict nonlinear effects of time on erosion rate, to quantitatively arrive 
at meaningful correlations in the laboratory, and to extrapolate to field 
prototypes. Figure 2 present theoretical prediction curves of relative 
intensity (or relative erosion rate) versus relative time for an attenuation 
exponent n = 2. The final equation used is also presented in figure 2. This 
theory needs four parameters to compute the erosion-rate-versus-time curve 
(table I). Long-term predictive efforts using the theory of erosion (ref. 10) 
were also attempted by Thiruvengadam (ref, 22). Thiruvengadam also modified 
his nomogram using the concept of erosion strength instead of strain energy 
(ref. 23). For design engineers, this nomogram may possibly be useful for the 
rough estimation of erosion rate under cavitation and liquid impingement 
erosion conditions. 


Hoff and Langbein Equation 

A simple exponential equation was proposed by Hoff and Langbein (ref. 11) 
incorporating the heterogeneous characteristics of impingement drops based on 
impact statistics. The proposed equation (fig. 3) requires only two parameters 
in order to compute the erosion rate as a function of time (table I). It 
should be noted, however, that by introducing a Poisson distribution into the 
method proposed by Heymann (ref. 8) or by introducing a distribution function 
into the original method proposed by Hoff and Langbein (ref. 11), the two 
methods are quite similar. 


Tichier and de Gee Model 

Tichler and de Gee (ref. 9), on the basis of the observation of two 
steady-state periods, have formulated an equation to predict the mean depth of 
erosion as a function of time. It was assumed that the erosion rate is 
relatively high and the surface is attacked uniformly during the first steady- 
state period. The surface is saturated with deep isolated craters and the 
erosion rate is relatively low during the second steady-state period. The 
final equation suggested for attenuation and the second steady-state period and 
definition of terms used are presented in figure 4. The equation needs five 
parameters to define the mean depth of erosion-rate-versus-time curve 
(table I). A graphical method to determine these parameters was presented by 
the investigators. 


Noskievic Formulation 

Noskievic (ref. 13) formulated a mathematical relaxation model for the 
dynamics of caviation damage of materials using a differential equation 
applied to forced oscillations with damping. This cavitation erosion model 
requires three parameters for the prediction of an erosion-versus-time curve 
(table I). Charts are presented for the use of this method (fig. 5), which 
simplifies the hurdle of going through lengthy equations and calculations. 

The experimental curve of relative cavitation damage (o) versus log cavitation 
exposure time (log t) has to be compared with curves in figure 5 by shifting 
in the direction of the log t axis until a curve of approximate match is 
found. This enables one to read out 6 and x, which results in a 3 value. 

Curve-Fit Approach 

Data for a large number of materials tested in both a rotating disk 
device and a magnetostriction oscillator have been analyzed in a manner that 
presents normalized cumulative average erosion rate versus normalized time 
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which brings the results to a universal curve fit (ref. 14). With a knowledge 
of four parameters (table I), it may be possible to correlate erosion data 
between the laboratory model and field devices. The agreement of the data 
analyzed from two previous investigations with entirely different experimental 
conditions not only showed similarities between cavitation and liquid— 
impingment erosion, but also reinforced the possibility of the unified nature 
of erosion. Correction factors for the incubation period and intensity of 
erosion are suggested. 


DISCUSSION 

Figure 6 presents normal ized-average-erosion-rate-versus-normalized-time 
curves for stainless steel tested in a rotating disk device (ref. 14). The 
data are normalized with respect to peak erosion rate and the time corre- 
sponding to this peak. Models proposed by Thiruvengadam (ref. 10), and Heymann 
(ref. 12), and a curve-fit normalization technique by Rao and Y'vjng (ref. 14), 
have also been presented on the curve. The methods proposed in references 10 
and 12 fit the data following the peak erosion rate. It is noted, however, 
that normalized time from 0 to 1 cannot be represented by any of the equations 
presented earlier, except the curve-fit approach developed at this laboratory 
(ref. 14). The methods proposed by Tichler and de Gee (ref. 9) and Noskievic 
(ref. 13) have not- been used, as the data considered for the analysis was not 
exposed too long and their plots represent a different dependent parameter. 

For individual materials, good results can be obtained at a single experimental 
condition with these two methods. It must be indicated, however, that many 
calculations are needed with these two methods. 

To check the general validity of the models and graphical approaches 
presented earlier, data reported for cavitation erosion (refs. 24 and 25) and 
liquid impingement (ref. 26 and 27) were analyzed. A typical set of plots is 
presented in figures 7 to 9 as normalized average erosion rate versus nor- 
malized time. It is evident that a material tested at a variety of conditions 
cannot be represented by a single method proposed earlier for long-term predic- 
tions. When individual groups of materials are considered, as in figure 9, the 
methods proposed by Thiruvengadam (ref. 10) and Heymann (ref. 12) are good. To 
show the involvement of calculations with equations proposed in reference 9, 
table II presents parameters necessary to calculate the mean-depth-of-erosion- 
versus-time curve for stainless steel, mild steel, and brass tested in a liquid 
impact device (ref. 26). In order to use this method, one must know the two 
steady-state periods. These may, however, not be available for most of the 
materials tested. The curve-fit approach suggested by this laboratory 
(ref. 14) produces a large scatter band to cover a wide variety of experi- 
mental conditions. This method, however, not only calculates erosion rates 
and times, but also the cumulative erosion. 

In addition to a curve-fit approach (ref. 14) suggested earlier to solve 
long-term predictions, a new characteristic law of average erosion rate versus 
cumulative erosion is also presented herein. 

Characteristic of Erosion-Rate-Versus-Cumulative-Erosion Curve 

Figure 10 presents a typical plot of cumulative average volume loss rate 
versus volume loss of mild steel tested in a rotating disk device. The experi- 
mental conditions are: pressure, 0.15 MPa (abs); velocity, 37.3 m/sec; 

diameter of the cavitation inducer, 25.4 mm; and diameter of the test specimen, 
63.5 mm. It appears that this curve has acceleration, peak rate, and decelera- 
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tion zones. The acceleration and deceleration zones may be represented by 
separate power-law relations. The equation for the acceleration zone is 
written as: 


V/t » AV n 

or V « ( At ) ^ ^ 1— n ) 

ORIGINAL PAGE & W 

OF POOR QUALITY 

where V « cumulative volume loss, mn^; t * exposure time corresponding to V, 
min; A * coefficient; and n « exponent. Differentiation of equation (2) with 
respect to t and after simplification results in 

dV/dt « V/ B l-n)t] 

(3) 

Similarly, the deceleration 

zone after the peak is represented as 

V/t « BV" m 

(4) 

or V = 

(5) 

where B *= coefficient, and m 
respect to t results in: 

* exponent. Differentiation of equation (5) with 

dV/dt = V/ Bi + m)tl 

VO/ 

= C t-P (1+m) 

(7) 

where C » B 1 /( 1+r1 ) and p = m/(l+m). The coefficients, exponents (slopes), 
and correlation coefficients obtained oy least-square fit are marked on 
figure 10. Equation (3) and (6) indicate that instantaneous erosion rate dV/dt 
for these two zones is a function of cumulative average erosion rate V/t. 
Further the ratios of these two rates are always constant. The intersection 
point for these two curves may be obtained by equating equations (1) and (4), 
i • e • , 

AV n « BV~ m 

(8) 

or V « (B/A) 1 /( n+m ) 

(9) 


The value of V in equation (9) corresponds to maxima on average-erosion-rate- 
versus-erosion curve, and values of (V/t) max and time corresponding to 
this peak may be obtained by using either equation (1) or (4). 

This study establishes that exponents n and m are almost equal, and 
that a power-law relation also exists between instantaneous erosion rate and 
exposure time. The advantage of this characteristic relation is that the 
values of ( dV/dt ) max and time corresponding to this peak may be calculated 
with only a few experimental points. It is generally observed that 
equation (1) terminates at (dV/dt ) max and deviates from the experimental 
points. To the knowledge of the present authors this type of power-law 
relationship has not been reported earlier. This relationship opens new 
avenues in erosion scaling and modeling efforts. 

CONCLUSIONS 

A brief overview of long-term cavitation erosion prediction equations and 
their capabilities is presented. Data analysis using cavitation and liquid 
impingement erosion data indicates that the normalized curve-fit approach 


suggested from this laboratory affords a better prediction for certain sets of 
data. For individual materials at one experimental condition, however, the 
methods proposed by Thiruvengadam and Heymann are good immediately following 
the peak erosion rate. 

A unique power-law relationship between average erosion rate and cumula- 
tive erosion is presented. It is believed that this relationship can solve 
some long-term modeling problems. 
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TABLE I. - PREDICTIVE MODELS, FORMULATIONS, AND PARAMETERS NECESSARY 
FOR COMPUTATION OF CAVITATION AND LIQUID-IMPINGEMENT 
EROSION VERSUS TIME CURVES 


Investigator 

Type of erosion 3 


Thiruvengadam 

CAV (nomogram) 

( 1 ) 

(refs. 19 
and 21) 


(2) 

Heymann 

CAV and LI 

U) 

{ref. 8) 

(elementary 
model ) 

(2) 

4 


LI (elaborated) 

(1) 


model ) 

(2) 



(3) 

Thiruvengadam 

LI and CAV 

(i) 

(ref. 10) 


(2) 



(3) 



(4) 

Hoff and 

LI (rain 

(1) 

Langbein 
(ref. 11) 

erosion) 

(20 

Heymann 
(ref. 12) 

CAV and LI 

(1) 

(2) 



(3) 


LI 

(1) 


( 4 ) 


Parameters needed for computation 

erosion intensity 
strain energy 


Nominal mean lifetime for original 
surface 

Standard deviation for original surface 
Nominal mean lifetime for substructure 
Standard deviation for substructure 

Delay time during which no failure occurs 
Mean of log-normal distribution on 
logarithmic time scale 
Standard deviation of log-normal 
attribution on logarithmic time scale 

Magnitude of instantaneous erosion rate 
at first peak, I^g* 

Time to attain first peak instantaneous 
erosion rate, t m 
Attenuation exponent, n 
Weibull shape parameter, a 

Maximum rate of erosion, e max 
Incubation period {intercept on time axis 
from straight line portion of erosion vs. 
time curve), t 

Mean depth of erosion at tangent point, 

yj 

Average erosion rate at tangent point, 

Rt 

Exponential constant, B 

Cumulative mean depth of erosion or 
material loss at tangent point, yj 
Normal component of impact velocity, V 0 
Volume of liquid impinging per unit area 
per unit time, U ? 

Generalized nondimensional erosion 
resistance parameter. N e 
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Investigator 


Parameters needed for computation 


Type of erosion 9 


Tichler and 

CAV 

(1) 

de Gee 


(2) 

(ref. 9) 


(3) 



(4) 




(5) 

Perelman and 
Denisov 

LI 

U) 

(ref. 17) 


(2) 



(4) 



(5) 



(6) 



(7) 

Thiruvengadam 

CAV and LI 

(1) 

(ref. 23) 

(nomogram) 

(2) 

Lichtarowicz 

LI (graph) 

(1) 

(ref. 20) 

(2) 

Noskievic 

CAV 

(1) 

(ref. 13) 


(2) 



(3) 

Rao and Young 

CAV and LI 

(1) 

(ref. 14) 


(2) 


Incubation time, t 0 
Resistance against cavitation erosion 
under hydrodynamic conditions, as occur 
in magnetostrictive oscillator R 5 
Mean depth of erosion at which effect of 
crater formation becomes manifest, r^ 
Proportionality constant, symbolizing 
increase in mean depth of erosion that 
would be necessary to form number of 
craters per unit area in final of craters 
per unit area in final steay-state 
period, < 

Ratio of rate of erosion in the final 
steady-state period to rate of erosion in 
first steady-state period, 

work done on microplastic deformations 
per cycle of load, due to the energy 
capacity of microvolumes of the material 
energy expended on fatigue fracture 
influence of the surface form 
kinetic energy of a stream of droplets 
energy in a steam of droplets absorbed 
during the incubation period 
specific energy of fracture determined 
from macroscopic fracture tests 
energy absorbed by the material during 
initial deformation 

erosion intensity 
erosion strengtu 

cumulative peak erosion rate 

time to reach cumulative peak erosion 

rate 

Cavitation property of material, a 
Cavitation strength of material or 
inner friction of material during 
plastic deformation, a 
Cavitation damage rate in developed 
period of cavitation attack, v s 

Peak cumulative average erosion rate 
Time to attain peak cumulative average 
erosion rate 
Incubation period 
Erosion resistance 


a CA V: cavitation erosion 

LI: liquid impingement erosion, cyl indrical /spherical drop or jet impact 

including jet with cavitation inducer. 
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TABLE II - PARAMETERS NECESSARY FOR TICHLER AND DE GEE MODEL 

‘iata source: 


Parameter 

Stainless Steel 

t 0> impacts 

252 x 10 3 

R c 

2.32 x 10- 3 

r b , m 

328 

R.o 

5.19 x 10- 4 

<?=R ro /R^ 

0.244 

a 

1.422 

6 

1.813 

r(t s ) 

369 


392 

k; 

75.62 


ref. 26 


Mild steel 

60/40 brass 

85 x 10 3 

189 x 10 3 

2,72 x 10~ 3 

4.20 x KT 3 

360 

620 

1.46 x 10“ 3 

3.44 x 10~ 4 

0.536 

0.082 


0.763 


3.395 

401 

722 

453 

849 
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OF POOR QUALITY 
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LIQUID IMPINGED PER UNIT 
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Figure 1. -Typical cumulatlvo erosion versus 
time curve, defining various terms usod 
and equations suggested, (Ret. 11). 
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Figure 2, - Theoretical prediction of the effoct of time 
on intensity of erosion when n * 2 (Ref. 9). 
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Figure 4. - Equation and definition of parameters used by rich ler and do Gee (ref. 81. 
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Figure 5. - Relative cavitation damage versus relative time curves 
(ref. 12). 
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(SEE FIG. 2) 
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Figure 6. - Normalized average erosion rate versus normalized time of stainless 
steel In a rotating disk device. 
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Figure 7. - Normalized average erosion rate versus normalized time of different materials - 
vibratory caviiation. 



Figure 8, - Normalized average erosion rate versus normalized time of L-605 alloy tested at different 
temperatures in liquid sodium. 



